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Introduction 

A  frequently  encountered  class  of  network  polymers  consists  of  polymer 
chains  crosslinked  via  a  small  monomer  crosslinkcr.  It  is  useful  to  be  able  to 
predict  network  parameters  such  as  gel  point  and  weight  fractions  of  soluble, 
pendant  and  elastically  effective  material.  Probability  models  can  be  developed  to 
accomplish  this,  subject  to  certain  simplifying  assumptions. 

The  purpose  of  this  paper  is  to  extend  previous  probability  models  of 
crosslinking  systems  to  include  polymer  chains  with  an  arbitrary  number  of 
different  types  of  monomer  repeat  units,  some  with  reactive  sites  and  others 
without.  A  Markov  model  of  chain  growth  is  assumed.  Certain  statistical 
properties  of  the  polymer  chains  are  needed  to  compute  network  parameters  of  the 
crosslinked  system.  These  properties  are  computed  using  the  classical  theory  of 
finite-state,  transient  Markov  chains;  it  is  not  necessary  to  resort  to  Monte  Carlo 
techniques  to  estimate  properties  of  such  general  chains. 

A  network  model  is  based  on  the  three  simplifying  assumptions  of  Flory  121 
and  Stockmayer  (81:  (i)  functional  sites  have  equal  reactivity,  (ii)  functional  sites 
react  independently  of  one  another,  and  (iii)  intramolecular  reactions  do  not  occur  in 
the  finite  species.  Using  statistical  properties  of  the  polymer  chains  and  the 
crosslinking  monomer,  the  network  model  computes  gel  point  and  network  parameters 
such  as  weight  fractions  soluble,  pendant,  and  elastically-effcctive  as  a  function  of 
conversion.  This  model  shows  the  effect  of  unequal  rates  of  homopropagation  and 
crosspropagation  in  the  chain-growth  reaction  on  these  network  parameters. 

A  Model  of  Polymer  Chain  Growth 

Polymer  chains  are  formed  from  several  reactants:  initiator,  1;  transfer 
agent,  TH;  monomers,  M,^,  k==l,  2.  .  .  .  m.  Chain  growth  occurs  in  tl’.rce  stages. 
During  initiation,  the  initiator  I  changes  into  two  activated  initiators.  I*.  The 
activated  initiator  then  reacts  with  an  unreacted  monomer,  Mj.,  resulting  in  a  short 
chain  consisting  of  a  reacted  initiator,  I^,  and  an  activated  monomer,  MJ., 
k  =  l,  2,  .  .  .  ,  m.  Chain  growth  consists  of  the  activated  chain-ends  reacting  with 
unreacted  monomers.  Chain  termination  can  occur  in  three  ways.  Termination  by 
combination  occurs  when  two  chains  with  activated  cnd-moncmers  combine  to  form 
one  completed  chain.  Termination  by  disproportionation  o>_curs  when  two  chains 


with  activated  end-monomers  terminate  but  do  not  combine,  resulting  in  two 
completed  chains.  Termination  by  transfer  occurs  when  a  live  chain-end  reacts  with 
a  transfer  agent,  TH;  the  hydrogen  atom  completes  the  chain  and  the  transfer  agent 
is  activated,  T*.  The  activated  transfer  agent  T*  then  initiates  a  new  chain  by 
reacting  with  an  unreacted  monomer.  This  chain  growth  reaction  can  be  represented 
as 


Initiation: 

I 

2  1* 

k‘1 

I*  + 

M. 

- 

l^M- 

k“ 

J 

T*  + 

T’-M* 

k^‘ 

J 

Growth; 

+ 

Mj 

- 

~m[m* 

Termination: 

~M*  + 

k^^ 

>J 

~M*  4- 

~Mj 

+  ~m;. 

TH 

- 

-MJ’H  -f  T*  , 

*^i 

for  i,  j  =  1,  2.  .  .  .  ,  m,  where  the  k's  represent  kinetic  rates.  This,  and  similar, 
chain  growth  reactions  are  described  by  Williams  (91,  Odian  (71,  Galbraith,  et  al.  [31, 
and  others.  A  system  of  kinetic  equations  based  on  the  above  scheme  with  second 
order  kinetics  can  be  defined  and  solved  (by  numerical  integration).  This  yields 
concentrations  (or  numbers)  of  1,  l'"  T*,  T*^,  H,  M|^,  Mj^,  k=l,  2,  .  .  .  ,  m, 

as  a  function  of  time.  This  system  consists  of  3m  4-7  differential  equations. 


It  is  assumed  that  the  above  system  is  stiff:  the  rate  at  which  activated 
initiator  is  introduced  into  the  system  is  much  less  than  the  rates  of  chain  growth 
and  termination.  This  means  that  the  chains  initiated  during  a  short  time  interval 
(of  length  At)  will  see  a  constant  environment  while  growing;  the  salient  aspects  of 
this  environment  arc  the  concentrations  of  Til,  Mj^  and  Mj.,  k-  1,  2,  .  .  .  ,  m.  This 
constant  environment  and  the  kinetic  growth  which  depends  on  just  the  last  unit  of 
a  live  chain  lead  to  a  Markov  model  of  chain  composition  for  the  chains  initiated 


during  the  interval.  Chains  from  other  intervals  will  also  have  Markov  models,  but 
with  different  parameters.  The  totality  of  chains  produced  will  correspond  to  a 
mixture  of  Markov  models.  Mixing  destroys  the  Markov  property;  thus,  the 
solution  approach  consists  of  computing  the  statistical  chain  parameters  needed  for 
network  analysis  for  a  succession  of  short  time  intervals  using  the  Markov  model 
and  parameters  for  that  interval,  and  then  taking  the  mixture  of  the  chains  and 
parameters.  The  number  of  chains  initiated  by  1*  and  T*  during  a  short  interval  of 
length  At  can  be  represented  by  Al*^  and  AT*^,  respectively.  These  values  and  the 
current  values  of  TH,  Mj^,  and  Mj^  uniquely  determine  the  polymer  chains  initiated 
during  the  time  interval. 


Instantaneous  Markov  Model  of  Chain  Composition 

Consider  the  polymer  chains  initiated  in  a  short  time  interval,  focusing  on 
the  portion  between  initiation  and  termination.  Each  position  can  be  filled  by 
different  elements.  The  initiating  elements  may  be  I  or  T.  Interior  elements  may  be 
Mj^,  k==l,  2,  ....  m.  The  terminating  element  may  be  H,  or  phantom  elements  D,  or 
Cj^,  k=l,  2,  .  .  .  ,  m,  where  D  represents  termination  by  disproportionation  and  Cj, 
represents  termination  by  combination  with  M|^,  k=l,  2,  .  .  .  ,  m.  Let 

X  =  (Xq.  Xj,  X.> . Xy)  represent  the  polymer  chain;  X  is  a  transient  Markov 

chain  with  initial  distribution  Tq  and  transition  probability  matrix  Pr(-  ,  ).  where 
the  subscript  “r”  represents  “right”,  i.e.,  the  chain  is  envisioned  growing  to  the 
right.  The  initial  distribution  is 


=  p(X(,  --  1)  =  Al’’/(AI‘’  +  AT*") 

tTqCT)  ■-=  P(Xq  -  T)  =  AT'’/(A1''  AT*"). 


The  transition  probabilities  PpCu,  v)  =  P(Xj^j  =  vlX^  =  u)  are  defined  for 
u  c  IM^,  .  .  .  ,  Mjni  I»  T)  and  v  c  IM^,  .  .  .  ,  Mm*  H,  D,  Cj,  .  .  .  ,  Cm)*  1  hey  are 
computed  from  kinetic  rates  and  the  instantaneous  concentrations  of  PH,  Mj^,  and 
MJ.,  k=l,  2,  .  .  .  ,  m.  For  example. 


PpH,  Mj) 


[M 


f 


This  model  describes  polymer  chain  structure  between  initiation  and  termination. 
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Now  consider  the  chains  which  terminate  by  combination.  A  chain  ending  in  a 
phantom  state  is  actually  combined  with  a  chain  component  ending  with  Mj^. 

This  component  is  viewed  in  the  “left”  direction,  i.e.,  the  structure  is  seen  looking 
toward  the  initiator,  rather  than  away  from  it.  The  Markov  property  still  holds 
but  the  index  set  for  the  chain  is  reversed.  Let  ?£(•  ,  •)  be  the  transition  matrix 
of  the  reversed  chain,  the  subscript  denoting  the  direction  “left,”  toward  the 
initiating  units.  The  reversed  transition  probabilities  P^lu,  v)  are  defined  for 
u  €  {Mj,  .  .  .  ,  Mn,,  H)  and  v  e  {M^  .  .  .  ,  Mnn  I,  T).  can  be  computed  from  t,o 
and  Pf.  Feller  [1]  and  Kemeny  and  Snell  (4]  are  recommended  references  for  the 
necessary  background  on  Markov  chains. 

The  initial  distribution  and  the  transition  matrices  Pj-  and  P£  could  be 
used  for  Monte  Carlo  generation  of  a  statistical  sample  of  chains  formed  during  the 
current  short  time  interval.  Note,  however,  that  chains  terminating  by  combination 
have  two  initiators  and  in  effect  will  be  counted  twice;  so  the  statistical  sample 
must  be  adjusted  accordingly. 

The  Network  Model 

The  chains  produced  in  the  above  reaction  are  crosslinked  by  reacting  them 
with  a  small  monomer  crosslinker,  Bg.  It  is  assumed  that  chain  monomer-unils 
k=l,  2,  .  .  .  ,  nf,  each  have  a  reactive  site,  denoted  as  an  A-site.  None  of  the 
other  chain  units  have  reactive  sites. 

Miller  and  Macosko  15]  present  a  network  model  for  crosslinking  polydispersed 
chains.  Their  equation  (107)  for  the  gel  point  is  valid  for  the  current  system: 

ac  =  (r(EsF^-l)  (EsFg-Ur^^^, 

where  r  is  the  stoichiometric  imbalance  and  expected  functionality 

of  a  chain  chosen  randomly  by  picking  a  reactive  site.  Their  analyses  is  based  on 
the  quantities  P(F^^^)  and  P(Fg^^)  which  are  solutions  of  the  equations: 

=  I  _  m  -f-  a  V  (P(FS“M)''  Ps(Pa-^='‘^ 
f-0 

=  1  -  r.x  -t-  ra  2  PgCFg  -  I -g). 

g-0 
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These  equations  are  also  valid  for  the  current  system. 

For  general  chains,  additional  modelling  is  required  for  computation  of  weight 
fraction  soluble,  pendant  and  elastically-effective.  From  any  repeat  unit  on  a  chain, 
it  is  important  to  know  the  network  structure  looking  “left"  and  “right”  along  the 
chain.  Let  Dg  represent  the  event  of  pendant  (dangling)  material  to  the  left.  Dp 
pendant  to  the  right,  and  -  Dg  Pi  Dp  pendant  in  both  directions.  For  a 
particular  repeat  unit  Mj,  the  probabilities  of  these  events  are 

P(n.iMi)  =  2  (PCF®*^^))*”  P(Fx=flM,) 

for  d  =  r,  £,  b;  where  Fr,  F^,  F^^  are  random  variables  representing  the  number  of 
functional  sites  seen  when  looking  “right,”  “left,”  and  in  “both”  directions  from  a 
repeat  unit  Mj,  i  =  1 ,  2,  .  .  .  ,  m;  and  from  1,  T,  and  H.  For  i  =1,  2,  .  .  .  ,  nf,  the 
monomer  units  have  A-sitcs;  thus 

P(Mj  soluble)  _  P(Fj“b  P(D^|M,) 

P(M,  pendant)  -  (1  -P(F°“^))  P(Dj,lM,) 

+  P(F5^^b(P(DrlM,)  +  P(D^,lMj)  -  2P(Djj,M,)) 

PCM-  effective)  =  1  —  P(Mj  soluble)  —  PCM^  pendant). 

For  i  —  nf -f- 1,  .  .  .  ,  m,  the  monomer  units  have  no  sites;  thus 

P(M,  soluble'  P(DjjlMj) 

P(M^  pendant)  =  P(DrlM,)  +  PCD^.M^)  -  2P  D^jlM,) 

P(Mj  effective)  1  —  P(M-  soluble)  —  P(M|  pendant). 

Similar  equations  hold  for  I,  and  II.  The  above  equations  can  be  combined  to 
compute  overall  weight  fraction  of  soluble,  pendant,  and  clastically-cffccti vc  chain 
material.  Using  the  above  equations  and  the  appropriate  theory  from  151,  network 
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weight  fractions  and  gel  point  can  be  computed,  provided  the  distributions  of 
functionality  to  the  right,  left,  and  both  directions  from  each  unit  are  available. 

Functionality  to  the  Right  and  Left 

Let  F£  equal  the  functionality  to  the  left  of  a  chain  unit,  M^,  say.  Consider 
the  Markov  chain  looking  left  which  ignores  Mj^  without  sites  (k=nf-fl,  •  •  •  .  m). 
Denote  the  chain  by  Y  =  (Yq,  Y^  .  .  .  ,  Y^);  where  Yq  equals  M-,  the  starting 
point,  and  Yy  equals  I  or  T,  the  terminations  looking  left.  The  transition  matrix 

r 

PgCu,  v)  is  defined  for  u  «  {M^,  .  .  .  ,  Mm.  H)  and  v  e  {Mj,.  .  .  ,  M^^p  I,  T},  and  can 
be  computed  from  using  the  concept  of  absorption  probabilities  for  Markov 
chains.  The  distribution  of  Fg  is  related  to  Y: 

P(F£  =  f!Mj)  =  P(Yf+|  =  1  or  T|Yq  =  M-), 

r 

which  can  be  computed  from  ?£.  Denote  this  conditional  distribution  as  JL(F£|Mj), 
for  i=l,  .  .  .  ,  m;  and  also  for  H. 

Now  consider  functionality  to  the  right,  Fp.  This  is  more  complicated 
because  of  the  possibility  of  termination  by  combination.  Let  F?  equal  the 
functionality  on  the  original  chain,  and  let  Fp  equal  the  functionality  on  the 
continuation  added  by  termination  by  combination.  Consider  the  Markov  chain 
Z  =  (Zq,  Zj,  .  .  .  ,  Zj-)  of  monomer  units  with  functional  sites  to  the  right  until 
termination;  the  termination  state  Z-p  may  be  H,  D,  or  Cj^,  k=-^l,  2,  ....  m.  The 

r 

transition  matrix  Pj-Cu,  v)  is  defined  for  u  t  (Mp  .  .  .  ,  Mm.  1.  T)  and 
V  e  {Mj,  .  .  .  ,  M^j.,  H,  D,  Cj,  ....  Cml  and  can  be  computed  from  Pp.  The 
conditional  distribution  of  F?  given  the  termination  state  is  related  to  Z: 

P(F?  =  flM,,  C^,)  P(Zf^.j  =  C^'iZq  =  MiVPfC^.Mj) 

where 

P(C^iMj)  =  2  =  ^k'^O  =  ^i^’ 

f*0 

r 

which  can  bo  computed  from  Pp.  Denote  this  conditional  distribution  as 
/.(F^iMj,  for  1,  k  ^--1,  2,  .  .  .  ,  m;  and  similarly  for  1  and  T,  initial  slates. 
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Suppose  the  additional  functionality  to  th'}  right  results  through  termination 
with  phantom  state  Cj^;  this  is  equivalent  to  looking  left  from  Mj^,  which  may  or 
may  not  contain  a  functional  site  itself.  So 

JLd+FglM,^)  k  =  l . nf 

JL(F?|Mj,  C^)  =  L(Fj\C^)  =  I 

I  l(FglM^)  k=nr+l,  ...,m 

Combining  the  functionality  of  the  two  parts  of  the  chain  to  the  right  gives 

m 

t(FrlM^)  =  '£  /.(FflMj.  C^)  i.(F?|C^)  P(C^!Mj) 

k=l 

-P  H)  P(HIMi)  +  KF^IMj.  D)  P(DiMj) 

where  represents  the  convolution  of  probability  distributions. 

For  chains  produced  in  a  short  time  interval,  F^  and  F^,  arc  conditionally 
independent  given  the  starting  segment.  Thus  the  sum  in  both  directions  is  obtained 
by  convolution: 


X(Fjj'M,)  =  X(FglMj)  JLlFrlM,). 

1  1  I  •  I  m* 

Mixtures  of  Polymer  Chains 

The  above  probability  laws  for  functionalities  to  the  right,  to  the  left  and  in 
both  directions  from  chain  repeat  units  are  for  the  chains  produced  in  a  short 
interval  of  time:  (t,  t-t  dt).  The  probability  law  for  the  totality  of  chains  is  the 
mixture  of  the  laws  for  the  short  intervals.  Suppose  /.^(FxlMj)  is  the  law  for 
chains  produced  in  (t,  t+dt)  and  dMj(t)  equals  the  number  of  monomers  reacted  in 
(t,  t+dt).  Then  ,  if  the  reaction  runs  over  (0,  S),  the  probability  laws  for 
functionality  are  given  by 

rS  rS 

KFjIM,)  =  l^CFjiMpdM.d)/  dM,(t) 

Jo  Jo 
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for  d=r,  £,  b,  and  i  =  l,  2,  .  .  .  ,  m.  Similar  equations  bold  for  1,  H,  T.  These  are 
the  distributions  needed  for  computation  of  weight  fractions  soluble,  pendant,  and 
effective.  The  distribution  of  1  needed  for  the  gel  point  is  obtained  by  mixing 
X,(Fjj|Mp,  i=l,  .  .  .  ,  nf. 


A  Computational  Example 


The  following  system  is  a  variation  of  one  in  15).  Polymer  chains  have  two 
types  of  repeat  units:  the  first  has  u  reactive  site,  the  second  does  not.  Elements 
have  the  following  masses:  =  146;  M2  =  142;  I  =  50;  T  ==  75;  H  ---  1 .  The 

kinetic  rates:  =  .01;  kj*  =  k^  =  2.0;  k|‘  =  k^'  =  4.0;  =  8.0; 


Iftc  _  .tC  _  Kie  _  _  f.  n.  |,>•u  _  i,«.u  _  •  lu  _  lu  _  a  n.  bii'  _ 

—  *^12  —  *^21  —  *^2”^  —  o.u,  —  *^12  —  *^21  —  *^22  —  *^21 

=  k®^  =  1.0.  (This  system  has  equal  rates  of  homopropagation  and 
crosspropagation.)  Moles  of  initial  reactants:  1  =  1.0;  TH  =  2.U;  —  50.0;  M", 

=  50.0.  The  reaction  went  for  a  time  duration  of  1.50  time  units.  The  solution  of 


,Br  _  ^gr  _  .  gr 


the  kinetic  equations  was  done  in  three  time  intervals,  using  a  simpic-minded 
difference  equation  aporoximation  with  a  step  size  of  0.001  lime  units.  The 
resulting  chains  had  a  number-average  number  of  7.43  monomer  units  and  a  number- 
average  functionality  of  3.72. 


These  polymer  chains  were  then  reacted  with  a  small  Bg  crosslmkcr.  This 
monomer  was  poiydispersed,  consisting  of  B„  b^,  and  B4  units  which  had  weights 
168,  478,  and  788,  respectively,  and  mole  fractions  0.3,  0.4  and  0.3,  respectively. 
The  mass  of  the  reactants:  mass  of  A-chains  =  1000.00,  mass  of  li-monomcrs  — 
600.00.  There  was  a  stoichiometric  imbalance:  r  =  .862.  The  site-average 
functionalities  are  Ev;(F^)  =  7,43  and  Es(Fj3)  =  3.20,  which  gives  a  gel  point:  a.^. 
=  0.286.  The  weight  fractions  of  soluble,  pendant,  and  elastically  effective 
material  are  shown  in  Figure  1  by  the  solid  lines  and  referred  to  as  System  1. 


In  order  to  investigate  the  effect  of  unequal  rates  of  homopropagation  and 
crosspropagntion,  two  additional  chains  systems  were  computed.  Both  systems  are 
identical  to  the  above  system  except  for  chain-growth  rates.  System  II  has  faster 
homopropagation:  k^^  -  ki^  =  1  •->  ^^2  ^  ~  faster 

crosspropagation:  k^[  =  k^^  =  0.1,  k^!^  =  k^[  =  1.9.  Both  of  these  systems 
yielded  chains  with  identical  number  average  length  (7.43)  and  functionality  (3.72)  as 
before.  In  System  11  (homopropagation  »  crosspropagation)  Eg(F^)  —  10.95 
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and  0.230,  a  sienificant  change  front  before.  In  System  111  (hotnopropagation 

C  crosspropagation)  Es(Fy^)  —  7.00  and  olq  =  0.296,  not  much  change.  The  weight 
fractions  for  System  II  are  shown  in  Figure  1.  The  weight  fractions  for  System  III 
do  not  differ  much  from  System  I  (horeopropagation  =  crosspropagation). 

Conclusions 

A  model  based  on  Markovian  chain  growth  has  been  developed  which  shows 
that  Monte  Carlo  methods  arc  not  required  for  analyzing  crosslinked  systems  of 
such  chains.  The  model  has  been  programmed  in  WATFIV  and  runs  quite  fast  (the 
examples  of  the  previous  section  each  required  17  seconds  to  compile  and  run  on  an 
IBM  4381).  Preliminary  investigations  show  that  unequal  rates  of  homopropagation 
and  crosspropagation  can  significantly  affect  network  parameters.  A  detailed 
description  of  an  extended  version  of  this  model  and  analysis  is  in  preparation  16). 

It  will  include  pre-gel  molecular  parameters,  multiple  types  of  unequally  reactive 
sites  on  chains,  additional  chain-forming  mechanisms,  and  other  topics. 
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Figure  1.  Weight  fractions  of  soluble,  pendant,  and  effective  material  for  two 
different  crosslinking  systems  as  a  function  of  conversion.  The  solid  lines 
represent  System  I,  an  ideal  system  (equal  rates  of  chain  homopropagation  and 
crosspropagation).  The  dotted  lines  represent  System  II,  a  system  with  much 
greater  homopropagation  rate  than  crosspropagation  rate.  The  lower  curves 
represent  weight  fraction  soluble.  The  differences  between  the  two  curves 
represent  weight  fraction  pendant.  The  diffcTcnccs  between  the  upper  axis  and  the 
upper  curves  represent  weight  fraction  of  elastically  effective  material. 
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